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ABSTRACT 



We present partial numerical results in the Minimal Supersymmetric Standard Model with 
soft breaking of supersymmetry, and radiative breaking of the electroweak symmetry. We impose 
the additional relation m& = m T at the GUT scale. For the special case of the strict no-scale model, 
in which global supersymmetry breaking arises solely from soft gaugino masses, we find that M t is 
expected to be lighter than ~ 128 GeV. A higher upper bound for M t of ~ 138 GeV is predicted, if 
. Mb is at the lower end of its experimental uncertainty. 

CO 

^ ' 1. Introduction 

Qh! In this paper, I report some preliminary results of work in progress, done in 

collaboration with E. Piard and D. Castaho. In the last few years, it has become 
apparent, using the ever increasing accuracy in the measurement of the strong cou- 
pling, that Supersymmetry (SUSY) affords an elegant means to achieve gauge cou- 
pling unification |^, [| at scales consistent with Grand Unified Theories (GUTs) 
H ^ [7], ^ . Whereas in the Standard Model (SM) the three gauge couplings unify 
"two by two" forming the "GUT triangle," in the simplest Minimal SUSY Extension 
of the SM (MSSM) these gauge couplings spectacularly unify at a point (within the 
experimental errors in their values). Given that the scale of unification in these mod- 
els is generally above the lower bound set by proton decay, these supersymmetric 
grand unified theories (SUSY-GUTs) have regained increasing interest. Constraints 
from the Yukawa sectors of such models have also yielded interesting predictions 
for various low energy parameters including the top quark mass || |10|, [ll], |12| . 

Crucial to these analyses has been the use of the Renormalization Group 
(RG) in extrapolating from the GUT scale to experimental scales. Here, we use of 
our previous RG analysis for the standard model[|13]], including the supersymmetric 



two loop [3 functions for the MSSM. Also soft symmetry breaking (SSB) terms are 
added. These lead to a nondegenerate superparticle spectrum as well as to the 
radiative breaking of the electroweak symmetry. A brief discussion of the effective 
one loop potential is presented. The boundary conditions at the unification scale in 



these minimal low energy supergravity models are discussed. Next the numerical 
procedure employed is described. The treatment of thresholds and the "special" 
form of the (3 functions needed is discussed next. Similar analyses have appeared in 
the literature @, [0| ^ |I7|. 



In the following, we only concern ourselves with the simplest schemes which 
have a chance of reproducing the data. If these are indeed those chosen by nature, 
they have some very distinct signatures. The top quark, although allowed by the 
standard model to be as heavy as 200 GeV, comes out much lighter. The lightest 
Higgs particle also comes out light, so light in fact that it may be very hard to detect. 
Another feature is that these models all predict the existence of a stable particle, 
the LSP, whose abundance can be calculated, based on the standard cosmological 
scenario. LSPs are typically too light to bring £1 to one, the preferred value of 
theorists. 

2. Minimal Supersymmetric Standard Model 

In the minimal supersymmetric extension of the standard model, every parti- 
cle has a supersymmetric partner, with spin differing by a half []18fl. Supersymmetry 



requires a second Higgs field with opposite hypercharge to the first as the super- 
potential cannot contain both a field and its complex conjugate. The second Higgs 
is also needed for both chiral U(l) and SU(2) global anomaly cancellation and to 
give its sector a mass. For renormalizable theories, the superpotential can have at 
most degree three interactions. The superpotential for the MSSM is (suppressing 
the SU(2) and Weyl metrics) 

W = uY u <b u Q + lYAdQ + #Y e $ d L + /x$ A . (1) 

where the hat indicates a chiral superfield and the overline denotes a left handed 
CP conjugate of a right handed field, V> = ia 2 ^* R - The usual Yukawa interactions 
are accompanied by new Yukawa interactions among the scalar quarks and leptons 
and the Higgsinos in the supersymmetric Lagrangian. There are also new gauge 
Yukawa interactions involving the gauginos. The new purely scalar interactions 
form the scalar potential which is positive definite in supersymmetric theories. The 
scalar potential will be discussed in a subsequent section. A remarkable aspect of 
supersymmetry is that all these new interactions require no new couplings. Table 
I displays the SU(3) x SU(2) x t/(l) quantum numbers of the chiral (all left handed) 
and vector superfields of the MSSM. 

The last term is put in by hand to avoid the exact Peccei-Quinn (PQ) sym- 
metry that the superpotential would otherwise exhibit. Although there are other 
ways to treat the PQ symmetry, including producing an axion, we leave the other 
possibilities to a future investigation. 

3. Minimal Low Energy Supergravity Model 



Since no super particles have been observed experimentally, supersymmetry, 
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if truly present in nature, must be broken. One way to accomplish this breaking is 
to add to the Lagrangian soft breaking terms. 

The most general soft symmetry breaking potential for the MSSM can be 
written (including gaugino mass terms) 
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+ E ( A « V?*i*uQi + ^Yfd& d Q 3 + A» YptedLs +h.c.) , (2) 
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where y sa «sino is the Majorana mass terms for the gaugino fields, A; (suppressing the 
group index), corresponding to U(l), SU(2), and SU(3), respectively. In the above, 
there are in principle sixty-three soft symmetry breaking parameters. This is hardly 
desirable to explain the standard model which already has eighteen parameters! 
It is clear that we need a further organizing principle. While this may be done 
by invoking extra symmetries such as flavor blindness (blandness), one attractive 
possibility is to couple the standard model to N = 1 supergravity (SUGRA). In 
the minimal low energy supergravity model, supersymmetry is explicitly broken by 
induced soft breaking terms of this form but with fewer parameters [|IT)| . 



We will assume that the spontaneous breaking of the local N=l supersym- 
metry is communicated to the "visible" sector by weak gravitational interactions 
from some "hidden" sector. This spontaneous symmetry breaking of supergravity 
manifests itself at low energy as explicit soft breaking terms of global supersymme- 
try. 

The SUGRA Lagrangian is characterized by two arbitrary functions of the 
fields, a real function (the Kahler potential) that determines the kinetic terms of the 
chiral superfields, and an analytic function, transforming as the symmetric product 
of the adjoint representation of the gauge group, that determines the kinetic terms 
of the gauge fields. In models with minimal kinetic terms for the chiral superfields, 
this leads to a common (gravitino) mass, m , for all the scalars of the model. The 
presence of non-minimal gauge kinetic terms implies non-zero masses, Mj, for the 
gauginos at the GUT scale, M x - By further assuming gauge coupling unification, 
we can take the three gaugino masses to be equal. Furthermore, the trilinear soft 



couplings A%, A l J, and are all equal to a common value A . With minimal chiral 
kinetic terms, the bilinear soft coupling B is related to A as B = A - m . This 
scenario has obvious, desirable features. First, it is very predictive since it has a 
few parameters accounting for thirty-one new masses. Second, the universal nature 
of the squark and slepton masses at M x helps to avoid the appearance of unwanted 
flavor changing neutral current (FCNC) effects. In fact, one could argue that the 
absence of FCNCs hints at a universal mass for the scalars. 

All of these couplings will evolve to different values under the renormalization 
group. The complete scalar potential now appears as 

V = V F + V D + V soft , (4) 

where Vf contains the potential contributions from the F-terms 

V F = \TlY u Q + ^d\ 2 + fiY d Q+iY e L + ^u\ 2 
+ \Y u Q<S> u \ 2 + \Y d Q<P d + Y e L<P d \ 2 

+ \%Y U $ U + lY d <P d \ 2 + |iY e $ d | 2 , (5) 
and Vd contains the potential contributions from the D-terms 

-,'2 



+ y ( Q f rQ + VtL + $tT$ u + ^ d r<f d f 

+ §■( Q^Q - S f A*S - d\*d f , (6) 
8 

where f = (n, r 2 , r 3 ) are the 517(2) Pauli matrices and A = (Ai, . . . , A 8 ) are the Gell-Mann 
matrices. In general, one must impose constraints on the parameters to avoid charge 
and color breaking minima in the scalar potential. Some necessary constraints have 
been formulated, such as 

Al < 3(m% + ml + m%J , (7) 
A% < 3(m 2 Q +m 2 d + m% 4 ) , (8) 
A\ < 3(m 2 L + m 2 e + m 2 ) . (9) 



However, these relations are in general neither sufficient nor indeed always 



necessary 



20]. Their derivation involves very specific assumptions about the spontaneous 



symmetry breaking. 



4. Radiative Electroweak Breaking 

An appealing feature of the models we are considering is that they can lead 
to the breaking of the electroweak symmetry radiatively ]2"T] , ^2] , ^ The one 
loop effective Higgs potential in these models can be expressed as the sum of the 
tree level potential plus a correction coming from the sum of all one loop diagrams 
with external lines having zero momenta 



Vx —loop 

(A) - V tree {A) + Afi(A) 



(10) 



The one loop correction is given by 



1 M 2 3 

AVl{A) = 64^ Str{A1(ln ^"2 )} 
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64 -^E(- 1 ) 2sp ( 2s P + 1 K( ln A^-9) - (") 



where M 2 is the field dependent squared mass matrix of the model and m p is the 
eigenvalue mass of the p th particle of spin s p . The tree level part of the potential is 

V tree (A) = m\{K)^ d + ml{K)^ u + ml{K){^ d + h.c.) 

+ ^( *£*,, ~ f + *1™ U + *$7*, f (12) 



where 

m?(A) = m| d (A)+ M 2 (A) , (13) 

mi(A) = m| ii (A) + / i 2 (A) , (14) 

ml(A) = B(AMA) . (15) 

The parameters of the potential are taken as running ones, that is, they vary with 
scale according to the renormalization group. The logarithmic term in the one loop 
correction is important in making Vi-i oop {A) independent of A to this order (up to 
non-field dependent terms). 

Given the low energy scale of electroweak breaking, we must use the renor- 
malization group to evolve the parameters of the potential to a convenient scale 
such as M z (where the experimental values of the gauge couplings are usually cited) 
thereby making this leading log approximation valid. The exact scale is not critical 
as long as it is in the electroweak range. If we define, 

-2 2 9AVl 

m * = m * + -Q^T > ( 16 ) 

with vi = Vd, v 2 = v u and 

dAVi 1 ^ ^_ tn , 1N 2/1 m 2 ^dml 



^^(-1)^(2^1)^^-1)5, 



then minimization of the potential yields the following two conditions among its 
parameters 

1 2 m\ ~ m\ t an2 P 
2 mz= tan 2 /3-l 

where m| = {g' 2 + g1)v 2 /2, v 2 = v 2 u + v 2 , and 



\ m z = T~^Tp i ' (18) 



B/j, = i(m 2 + 7™?,) sin 2/3 , (19) 



where tan/3 = v u /v d - 



Although results based on the tree level potential cannot always be trusted, 
one can still use it to get some idea under what conditions electroweak breaking 
occurs. The renormalization group evolution of m| can be such that it turns 
negative at low energies, if the top quark mass is large enough, whereas m| runs 
positive. From Eq. (TT^), the scale at which breaking occurs is set by the condition 

m?(A 6 ) m 2 2 {A b ) - mj(A h ) = . (20) 

If the free parameters are adjusted properly, then the correct value of the Z n mass 
(M z = 91.17 GeV) can be achieved. 

In the tree level analysis, there is another critical scale that must be consid- 
ered. It is evident from Eq. (|T^ ) that the potential becomes unbounded from below 
along the equal field (neutral component) direction, if 

m?(A s ) + m 2 2 (A s ) < 2m^(A s ) . (21) 



Since m\m\ -m\>§ implies m\ + m\ > 2ml, condition can only occur at scales 
lower than condition (pCf), so A s < A b . From this analysis, one concludes that the tree 
level vacuum expectation values (VEVs) of the scalar fields obtained by minimizing 
the potential are zero above A b , and grow to infinity as one approaches A s where 
the potential becomes unbounded from below. It follows that the appropriate scale 
at which to minimize the tree level potential and evaluate the VEVs is critical. 
This scale must be such that the one loop corrections of the effective potential 
may safely be neglected. Only at such a scale can the tree level results be trusted. 
However, there is more than one scale involved, and therefore, it is difficult if not 
impossible to find a scale at which all logarithms may be neglected. Indeed, the use 
of tree level minimization conditions to compute the VEVs at an arbitrary scale 
(e.g., A = Mz) leads to incorrect conclusions about the regions of parameter space 



that yield consistent electroweak breaking scenarios ||25|| . When AVi is included, 
however, the value of A is not critical as long as it is in the neighborhood of M z - 

Ref. ]25|] gives a prescription for arriving at a scale (A) at which the tree 



level and the one loop effective potential results for the VEVs agree. Three qualita- 



tively different cases are considered. Following Ref. P5fl , we let M S usy parametrize 
the superparticle thresholds, then the cases can be characterized by the orderings: 
(a) Msusy < A s < A 6 , (b) A s < M SUSY < A b , and (c) A s < A b < M S usy- In each case, 
the prescription is to take A = m&x{M S usY, A s }. Two cases deserve special mention. 
Case (a) cannot be handled using the tree level analysis because v u ,v d — » oo near A s . 
Fortunately, phenomenological bounds rule this case out anyway. In case (c), there 
is actually no electroweak breaking. For scales below Msusy, the superparticles have 
decoupled and the effective theory is not supersymmetric. Therefore, the running 
mass parameters of the potential freeze into their values at A = M S usy at which scale 
there is no electroweak breaking. Finally, it must be emphasized that the apparent 
violent behavior of the VEVs with scale in the tree level analysis is an artifact of 
the approximation. The only physical potential is the full effective potential, and 
it either breaks electroweak symmetry or not. If it does, then the scalar fields have 



non-zero VEVs, and these VEVs are non-zero over all scales varying according to 
the anomalous dimension of their respective scalar fields (in the Landau gauge) . 

In this paper, we do not rely on the tree level analysis, rather we incorporate 
the one loop corrections. We include the dominant contributions from the third 
generation, that is, those of the top and stop, bottom and sbottom, and tau and 
stau pB|, B7[ . The one loop effective potential is constant against the renormalization 



group to this order around the electroweak scale. We choose the Z° mass as the 
scale at which to evaluate the minimization conditions. Eqs. (|T8| ) and (|T9| ) can be 
written 

^ {Mz) ~ tan 2 /?-l 2 (22) 

m? + m 2 . ) sin 23 

B ^ = ' (23) 

where m% ud = m% n + dAVi/dvl d and used to solve for n{M z ) and B(M Z ) given the 
value of all the relevant parameters at M z . We note that the form of Eq. ( |22| ) does 
not fix the sign of (i, and a choice for its sign must be made (/i is multiplicatively 
renormalized) . The right hand sides of these equations implicitly involve the VEV 
at M z - In a consistent scenario it would have the value of v(M z ) = 174.1 GeV. 
If the parameters are such that \i 2 < 0, then the scenario is inconsistent and the 
electroweak symmetry fails to be broken. 

5. Boundary Conditions at M x 

In this paper, as in the previous one |13] we work in the modified minimal 
subtraction scheme (MS) of renormalization. The parameters of the Lagrangian are 
not in general equal to any corresponding physical constant. For example, in the 
case of masses, except for those of the bottom and top quark (see ]T3|]), all other 
physical masses, M, will be determined from their corresponding running masses by 
the relation 

M = m(A)| A=M . (24) 

This equation is easily solved in the course of an integration of the RGEs for the 
different masses by noting the scale at which it is valid. 

Because SUGRA models make simplifying predictions about the soft pa- 
rameters at the unification scale, we initiate the evolution of the renormalization 
group equations at this scale. It has been demonstrated that the introduction of 
supersymmetry leads to gauge coupling unification at approximately ~ 10 16 GeV. 
Therefore we take M x = 10 16 GeV, and evolve down to 1 GeV, the conventional scale 



at which the running quark masses are given |28|] . 



At the unification scale, M x , all the scalars will have a common mass 



m Q ,{M x ) 



= m Ui {M x ) = m di (M x ) = m Li (M x ) 

= m ei {M x ) = m^ u (M x ) = m^ d (M x ) = m , 



(25) 



Table 2: Models 





A 


m 


mi/ 2 


So 


Strict JNo-scale 








free 





No-scale 








free 


free 


String Inspired 





free 


free 






as will the gauginos 

Mi (M x ) = M 2 [M x ) = M 3 (M x ) = mi/a ■ (26) 
The prefactors of the trilinear soft scalar terms are equal at M x 

A* (M X ) = A'i (M x ) = A« (M X ) = A . (27) 

Also we define the bilinear soft scalar coupling and the mixing mass at M x by 
B(M X ) = B , and ft(M x ) = Mo- 

Furthermore, to constrain the parameter space, we will take the bottom and 
tau masses equal at M x 

m b (M x ) = m T {M x ) . (28) 

This being the best motivated mass relation in supersymmetric grand unified the- 
ories 



29 



The present analysis restricts itself to three subclasses of soft symmetry 
breaking models. These have various soft parameters equal to zero at M x . The 
first class of models follows from the no-scale model [[J0| and has A = ™ = B = 
(strict no-scale model). In these models, only gaugino masses provide global 
supersymmetry breaking. The second class is the less constraining no-scale case 
that has only A = m = 0. A third class we consider with A = B = comes from 
string derived models. Table II lists these three possibilities. In both strict no-scale 
and string inspired cases, we must have B = 0. However, since for us B is an output 
rather than an input variable, B = results must be inferred from its behavior upon 
varying other input parameters. 

6. Numerical Procedure 



In this work as in Ref. []T3[], we will use routines, based on the "shooting" 
method, to solve systems of nonlinear equations. The method involves making a 
guess for the solution, then assessing its merits based on how well the equations 
are satisfied, given some tolerance. The process is optimized and iterated until the 
routine converges on a solution. 

As discussed previously, we start out runs at M x at which scale we can make 
simplifying assumptions about the soft breaking terms based on various SUGRA 
models. This requires that we use the solution routines to consistently find the M x 
values of all known low energy parameters such as lepton and quark masses and 



mixing angles and gauge couplings. This amounts to solving for sixteen unknowns 
(nine masses, three angles and a phase, and three gauge couplings). Alternatively, 
we could start our runs at M z or 1 GeV; however, this now requires solving for 
sixty-three unknowns (the values of the soft breaking terms at low energy) that 
must evolve to just four different values at M x . The efficiency of the former method 
is obvious. 

There are seven free parameters in the model we consider. These are A , B , 
m o, mx/2, Ho, tan/3, and m t . The two minimization constraints (|18|) and ( |19D reduce 
this set to five, which are taken to be A , mo, m 1 / 2 , tan/3, and m t . In the present 
framework, B and ho will be determined using the numerical solutions routines in 
conjunction with the minimization of the one loop effective potential at M z in the 
process of evolving from M x to 1 GeV. Minimization at M z will give B(M Z ) and 
n(M z ). To arrive at B and ho (their corresponding values at M x ), we employ the 
solution routine as follows. A guess for B and ho is made at M x and then the 
parameters of the model are run to M z at which scale the evolved value of B is 
compared to the minimization output value for B at M z . The same is done for h- 
If the compared values agree to some set accuracy, then B and ho are the required 
values. Other analyses that also extract B(M Z ) and ii{M z ) simply evolve these two 
parameters via their renormalization group equations back to M x to find B and 
Ho relying on their near decoupling from the full set of RGEs. We note that the 
sign of h is n °t determined from the minimization procedure, thus we must make a 
choice for it. To constrain the parameter space further, the bottom quark and tau 
lepton masses will be taken equal at M x . This equality is a characteristic of many 
SUSY-GUTs. This constrains the model to four free parameters, A , mo, my 2 , 
and tan/3. Demanding that m b (M x ) = m T (M x ) and achieving the correct physical 
masses for the bottom quark and tau lepton fixes the mass of the top quark which 
affects the evolution of the bottom Yukawa significantly. We shall assume gauge 
coupling unification, an assumption which appears reasonable when one considers 
SUSY models with SUSY breaking scales & 10 TeV. 

In a complete treatment, the solution routines would be used to find the 
precise (similar) values of a\, a 2 , and a 3 at M x that will evolve to the experimentally 
known values at M z , however this increases the CPU time considerably. We shall 
therefore sacrifice some precision in their M z values by taking them exactly equal 
at M X - This is already a theoretical oversimplification since one does not expect 
the gauge couplings to be exactly equal due to threshold effects at the GUT scale. 
We find that for all cases we have studied, the common value aj" 1 (Mj) = a 2 1 (M x ) = 
CX3 1 (M X ) = 25.31 leads to errors no bigger than 1%, 5%, and 10% in ai(M z ), a 2 {M z ), and 
olz{M z ), respectively. This is not so bad considering that the (combined experimental 
and theoretical) errors on a 3 (M z ) from some processes can be as large as 10% (13| . 

It is well known that there is a fine tuning problem inherent in the radiatively 
induced electroweak models. For certain values of the parameters, the top quark 
mass must be tuned to an "unnaturally" high degree of accuracy to achieve the 
correct value of M z . This problem is generally handled by rejecting models that re- 
quire "too much" tuning. The amount of tuning is usually defined quite arbitrarily. 



The usual procedure is to define fine tuning parameters 



9MI . 

z 1 (29) 



^ 'Ml dx\ 1 ' 

where Xi are parameters of the theory such as m , m\j2, Mi or m t . One then demands 
that the e* be less than some chosen value that is typically taken to be 10. 

We have analyzed to some extent the differences in using the tree level vs. 
one loop effective potential. The basis for the "theoretical" fine tuning problem can 
be seen, if one makes some simplifying assumptions, in the dependence of M w on 
the top quark Yukawa coupling, y t , [|31 



M w ~ M x e- 1/y < . (30) 

We remark that this fine tuning problem is exacerbated if one uses only the tree 
level analysis of the potential. The vacuum expectation value coming from the 
minimization conditions of the tree level potential changes rapidly from to infinity 
over the interval (A s , A b ). Using the prescription of Ref. [p5| for the scale A at which to 
adequately minimize the tree level potential, to extract v(A), and thereby to arrive 
at a value for Mz, one finds that although a small variation in y t {M x ) may lead 
to a small variation in A b , the steepness in the tree level VEV can lead to a large 
variation in the value of v(A) and therefore in M z . Hence, in the tree level analysis, 
solutions which may be within the bounds of the "theoretical" fine tuning may 
nevertheless display a fine tuning aspect because of this "tree level" fine tuning of 
A. However, our use of the one loop effective potential (i.e., including AVi) stabilizes 
the VEV around the M z scale and this particular fine tuning goes away. The true 
VEVs depend on scale through wave function renormalization effects which are 
never large as can be seen from the form of the renormalization group equations for 
the VEVs. 

In this analysis, we shall also reject solutions based on fine tuning consid- 
erations; however, our method differs somewhat from the usual one in that it is 
incorporated in the solution routine described above. The routine is an iterative 
one which determines the convergence properties of the solution. Very slow conver- 
gence reflects an inherent fine tuning. Therefore, if the convergence is too slow, we 
will reject the solution. Effectively we are rejecting any solution which the computer 
cannot pinpoint within an allotted number of iterations. 

Given values for A , mo, m 1 / 2 , tan/3, and sign(/z), the solution routines search 
for the values of v(M x ), mb, T (M x ), m t (M x ), B , and fi . The process by which B and 
fi are found was described above. The remaining three parameters are determined 
similarly. The routine makes a guess for v(M x ), m b ^(M x ), and m t (M x ), then the 
full renormalization group equations are evolved to 1 GeV calculating superparticle 
threshold masses in the process and minimizing the one loop effective potential 
at M z - The merits of the guess for v(M x ), m b ^(M x ), and m t (M x ) is assessed by 
comparing the resulting values of Mz, m T (l GeV), and m b (l GeV) with the expected 
ones. The process is iterated until the correct values are achieved to within some 



tolerance. 



7. Thresholds 

In the minimal low energy supergravity model being considered, the super 
particle spectrum is no longer degenerate as in the simple global supersymmetry 
model in which all the super particles are given a common mass, Msusy- In the 
simple case, one makes one course correction in the renormalization group evolution 
at Msusy- In the model with soft symmetry breaking, the nondegenerate spectrum 
should lead to various course corrections at the super particle mass thresholds. To 
this end, the renormalization group functions must be cast in a new form which 
makes the implementation of the thresholds effects (albeit naive) evident. Since 
the MS renormalization group equations are mass independent, particle thresholds 
must be handled using the decoupling theorem fl52| , and each super particle mass 
has associated with it a boundary between two effective theories. Above a particular 
mass threshold the associated particle is present in the effective theory, below the 
threshold the particle is absent. 

The simplest way to incorporate this is to (naively) treat the thresholds as 
steps in the particle content of the renormalization group (3 functions. This method 
is not always entirely adequate. For example, in the case of the SU{2) gauge coupling 
there will be scales in the integration process at which there are effectively a half 
integer number of doublets using this method. We believe, nevertheless, that this 
method does yield the correct, general behavior of the evolution. It is a simple 
means of implementing the smearing effects of the non-degenerate super particle 
spectrum. The determination of the spectrum of masses is done without iteration 
as is common in other analyses. Our method deduces the physical masses by solving 
the equation m(A) = A for each superparticle in the process of evolving from M x to 1 
GeV. The usual iterative method requires several runs to find a consistent solution. 

8. Analysis 

The tremendous computing task involved in analyzing the full parameter 
space of the soft symmetry breaking models, using the methods described as de- 
signed, would be far too time consuming given the computing facilities available 
to us. Therefore, in the following analysis, some simplifications will be made in 
the procedural method. First, only the heaviest family of quarks and leptons will 
have non-zero mass. This will cut down on the CPU time required for the solution 
routines to consistently find their values at M x - Second, as stated previously, the 
value of the strong coupling at M z will be allowed to vary from its central value 
of .113 by at most 10%. This translates into a similar error in the bottom quark 
mass. Third, the allotted number of Runge-Kutta steps, involved in numerically 
integrating the renormalization group equations, will be cut down to ~ 100. 

Our method involves four input parameters A , m , mi/ 2 , tan/3 (and the 
sign(^)). The output is i? , z^o, M t , M h , and all the masses of the extra particles 



associated with the MSSM. Efficient use of CPU time required that we proceed 
as foilows. For a given model, our initial exploration of the parameter space was 
performed in a coarse grained fashion, tan/3 was most commonly coarse grained as 
2, 5, and 10, with only some runs involving higher values (e.g., 15, 20). The other 
three input parameters were varied in steps of 50 and 100 GeV. Values larger than 
~ 500 GeV were rarely ever used. We subsequently narrowed down on the allowed 
hyperenvelope by fine graining around the edges of the expected region (based on 
the coarse graining results). 

Our raw data consists of those runs which satisfy the following two crite- 
ria. The first is consistent electroweak breaking; that is, the correct value of M z is 
achieved from the minimization of the one loop effective potential with fi 2 (Mz) > 0. 
The second criterion is no fine tuning, as implemented in our method (see Sec- 
tion 6). The solution routines employed return a numbered code representing the 
convergence properties of the solution which we use to screen the runs. 

The raw data is then progressively filtered based on three physical con- 
straints. R-parity is a discrete symmetry which distinguishes particles from super- 
particles by assigning a +1 to particles and a -1 to the superparticles. Conserved 
R-parity requires the existence of a stable lightest supersymmetric particle (LSP). 
Astrophysical considerations indicate that the LSP must be neutral and colorless. 
Cosmological considerations based on the LSPs contribution to the density of the 
universe indicate that it must have a mass less than ~ 200 GeV |3^, [53], |35|]. First, 
points in parameter space that lead to LSPs other than neutralinos with masses less 
than — 200 GeV are cut. Second, flavor changing neutral current bounds are used 
to reject runs in which intergenerational splitting of squark and slepton masses is 
too large. Third, experimental limits on the masses of the superparticles are used 
as another criterion to reject runs. 

9. Results 

In the following, we discuss the results of our analysis of the three classes of 
models listed in Table II. Because we have performed a coarse grained study, our 
results should only be considered qualitatively valid. Based on these, we hope to be 
able to ascertain the general trends in the data, and make some general predictions 
about the feasibility of the models considered. 

Because the GUT inspired constraint Eq. (^8p is enforced in this analysis, 
the results will depend on the mass of the bottom quark. Namely, lower bottom 
quark masses require larger values of the top quark mass to satisfy this relation. 
Most results will be reported for the case m fe (l GeV) = 6.00 GeV, but lower mass (5.70 
GeV) and higher mass (6.33 GeV) cases were also studied. The running value of 6.00 
GeV for m b (l GeV) corresponds to a physical bottom mass M b = 4.85 ± .15 GeV, with 
the uncertainty coming from the error in the strong coupling, as discussed above. 

9.1. No-scale Case 

In all cases considered, the mass of the LSP, when it is a neutralino, is 



observed to be correlated with the value of m l / 2 . Therefore, we find that my 2 
cannot be taken too large (< 400 GeV). In the no-scale case, we present plots of 
M t vs. m 1 / 2 for three values of m b (l GeV) and containing all points satisfying the 
various criteria outlined in Section 8. Figure 1 is such a plot for m b (l GeV) = 6.00 
GeV. The "right edge" of the envelope is defined by points whose neutralino LSPs 
are just slightly heavier than the lightest charged superparticle (usually a f R for us). 
The "left edge" defines the threshold of consistent electroweak breaking. The top 
and bottom edges are set by the requirement that Eq. (|28| ) hold. The range of tan/3 
considered leads to definite lower and upper bounds on M t || |10|, [11], [l^ . 



In Figs. 2 and 3, we display similar plots with lower (m b (l GeV) = 5.70 GeV) 
and higher (m b (l GeV) = 6.33 GeV) bottom quark masses. The previously noted 
dependence of the data on M b is evident from these figures. In Fig. 3, we change 
the value of the bottom quark mass to illustrate the dependence of our results on 
its value. 

From these, we note that the available range of mi/ 2 decreases with increasing 
M b . Fig. 1 indicates that 190 GeV ^ m l/2 & 265 GeV. We can draw no conclusions, 
however, about the value of tan/3. We can conclude from the m b (l GeV) = 6.00 GeV 
case that 85 ^ M t < 132 GeV. From similar plots involving M h , we conclude that 
25 ^ M h < 78 GeV. In Fig. 4, we display the familiar dependence of M t on tan/3 for a 
particular value of mi/ 2 in the allowed range. 

The official experimental lower bound on the top quark mass is 108 GeV. 
The figures indicate that the top quark cannot have a mass greater than ~ 132 GeV 
in this model, if m&(l GeV) = 6.00 GeV. This upper bound is raised to ~ 160 GeV if 
m b (l GeV) = 5.70 GeV, and the model is ruled out, if m b (l GeV) = 6.33 GeV. Figure 
5 is similar to Fig. 1, but we have chosen sign of /i negative. The allowed region is 
displaced down with respect to the positive fi case with the upper bound on M t now 
113 GeV and very close to the experimental limit. 

9.2. Strict No-scale Case 

The results of the strict no-scale case must be interpolated from the no-scale 
results, because B Q is not an input parameter in our procedural method. Therefore, 
we plot in Fig. 6 M t vs. B and deduce the M t bounds from slicing the data along 
B = 0. Admittedly, the relatively small number of points makes the perimeter of 
the region unclear in some areas. We get approximately 90 ^ M t ^ 127 GeV. For 
M hl we find 40 ^ M h <, 78 GeV with an uncertainty in the lower bound of ~ 10 GeV 
due to lack of definition in the lower end of the envelope. Inspection of the data 
indicates that 3.5 <, tan/3 < 9. Thus, it appears that tan/3 cannot be too small or too 
large to accommodate the strict no-scale case. 

Finally, in Table 3, we display the spectrum of superparticle masses for a 
representative strict no-scale scenario with tan/3 = 8.3 and m 1 / 2 = 240 GeV. The 
masses are in GeVs. There are several interesting features. For instance, the LSP 
is the photino in this case and has a mass of 92 GeV. It is too light to account for 
the dark matter. The top quark mass is just above the experimental limit at 126 
GeV, and the Higgs boson mass is 77 GeV, which underlines one of the signatures of 



Table 3: Strict No-scale Scenario 
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these scenarios: the Higgs is typically light, so that it will be found in two photon 
modes at SSC. At any rate, it illustrates the predictive power of this possibility. 

9.3. String Inspired Case 

As in the strict no-scale, our results for the string inspired case necessitates 
interpolating A = data to B Q = 0. Once again we find the perimeter of the allowed 
region is not well defined everywhere. Hence, our results are only qualitative. Figure 
7 is similar to Fig. 6, but in this case we only fix A = 0. The strict no-scale case is a 
special case of the string inspired one, therefore the 127 GeV top quark mass upper 
bound is not expected to decrease but rather to increase in this case. Slicing along 
B Q = yields 85 <, M t <, 140 GeV. Similarly, for the light Higgs we get 40 <, M h <, 80 
GeV. The data indicates in this case, as in the strict no-scale case, that there is a 
lower bound on tan/3 of — 3. 

10. Conclusions 

Minimal low energy supergravity models were considered. It is quite remark- 
able that for as simple a breaking of supersymmetry as that offered by the strict 
no-scale model, we can reproduce the standard model results, including the appeal- 
ing feature of automatic breaking of the electro-weak symmetry. The importance of 
the value of the top quark mass to these schemes cannot be understated, since we 
found in our study of specific models, rather restrictive upper bounds for the top 
quark mass. No-scale models in which only gaugino masses provide global super- 
symmetry breaking yield top quarks with masses less than ~ 127 GeV. The results 



are sensitive to the value of the bottom quark mass. Lower bottom quark masses, 
within the experimental uncertainty, lead to higher top quark upper bounds. In 
these models, the ratio of vacuum expectation values of the two Higgs fields is 
expected to be larger than ~ 70°. 

Although the perimeter of the allowed regions were often fuzzy, we could, 
nevertheless, draw some general conclusions from our results. For all our runs, with 
no restrictions on the soft terms, we find for the top quark M t < 186 GeV and for 
the light Higgs boson M h <, 100 GeV . 
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12. Figure Captions 

Figure 1. Plot of M t vs. m 1/2 for tan/? = 2,5, 10, 15 and m b (l GeV) = 6.00 GeV. 

Figure 2. Same as Fig. 1 but with m b (l GeV) = 5.70 GeV. 

Figure 3. Same as Fig. 1 but with m b (l GeV) = 6.33 GeV. 

Figure 4. Plot of M t vs. tan/3 in the no-scale case with m 1 / 2 = 240 GeV and 
m b (l GeV) = 6.00 GeV. 

Figure 5. Same as Fig. 1 but with sign(/i) = -. 

Figure 6. Plot of M t vs. B for points with A = m = used to interpolate to 
B = and deduce the M t range in this case. 



Figure 7. Same as Fig. 6 but for points with A Q = only. 



